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We propose a new theory to explain the formation of spiral arms and of all 
types of outer rings in barred galaxies. We have extended and applied the tech- 
nique used in celestial mechanics to compute transfer orbits. Thus, our theory 
£Nj ' is based on the chaotic orbital motion driven by the invariant manifolds as- 

sociated to the periodic orbits around the hyperbolic equilibrium points. In 
, particular, spiral arms and outer rings are related to the presence of hete- 

C ~ ~) ■ roclinic or homoclinic orbits. Thus, R± rings are associated to the presence 

of heteroclinic orbits, while R1R2 rings are associated to the presence of ho- 
moclinic orbits. Spiral arms and R2 rings, however, appear when there exist 
, neither heteroclinic nor homoclinic orbits. We examine the parameter space 

% ^ 1 of three realistic, yet simple, barred galaxy models and discuss the formation 

■ of the different morphologies according to the properties of the galaxy model. 

The different morphologies arise from differences in the dynamical parameters 
of the galaxy. 

Keywords: galactic dynamics - invariant manifolds - spiral structure - ring 
structure 



1. Introduction 

Bars are very common features in disk galaxies. According to Eskridge et 
al. [1] in the near infrared 56% of the galaxies are strongly barred and 6% 
are weakly barred. A large fraction of barred galaxies show either spiral 
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arms emanating from the ends of the bar or spirals that end up forming 
outer rings (Elmegreen & Elmegreen [2]; Sandage & Bedke [3]). 

Spiral arms are believed to be density waves (Lindblad [4]). Toomre [5], 
finds that the spiral arms are density waves that propagate outwards to- 
wards the principal Lindblad resonances, where they damp. So other mech- 
anisms for replenishment are needed (sec for example Lindblad [6] ; Toomre 
[5,9]; Toomre & Toomre [7]; Sanders & Huntley [8]; Athanassoula [10] for 
more details). Rings have been studied by Schwarz [11-13]. The author 
studies the response of a gaseous disk galaxy to a bar-like perturbation. 
He relates the rings with the position of the principal Lindblad resonances. 
There are different types of outer rings and they can be classified accord- 
ing to the relative orientation of the principal axes of the inner and outer 
rings (Buta [14]). If the two axes are perpendicular, the outer ring has an 
eight-shape and it is called Ri ring. If they are parallel, it is called R 2 ring. 
There are galaxies where both types of rings are present, in which case the 
outer ring is simply called R1R2 ring. 

Our approach is from the dynamical systems point of view. We first 
note that both spiral arms and (inner and outer) rings emanate from, or 
are linked to, the ends of the bar, where the unstable equilibrium points of 
a rotating system are located. We also note that, so far, no common theory 
for the formation of both features has been presented. We therefore study 
in detail the neighbourhood of the unstable points and we find that spiral 
arms and rings are flux tubes driven by the invariant manifolds associated 
to the periodic orbits around the unstable equilibrium points. 

This paper is organised as follows. In Sec. 2, we give the characteristics 
of each component of the model and the potential used to describe it. In 
Sec. 3, we give the equations of motion and we study the neighbourhood of 
the equilibrium points. In particular, we give definitions of the Lyapunov 
periodic orbits, the invariant manifolds associated to them, and of the ho- 
moclinic and heteroclinic orbits. In Sec. 4, we present our results and in 
Sec. 5, we briefly summarise. 



2. Description of the model 

We use a model introduced in Athanassoula [15] that consists of the su- 
perposition of an axisymmetric and a bar-like component. The axisymmct- 
ric component is the superposition of a disc and a spheroid. The disc is 
modelled as a Kuzmin- Toomre disc (Kuzmin [16]; Toomre [17]) of surface 
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density S(r) (see also left panel of Fig. 1): 

V 2 / 2\-3/ 2 

where the parameters Vd and set the scales of the velocities and radii 
of the disc, respectively. The spheroid is modelled using a spherical density 
distribution, p(r) (Eq. 2), characteristic for spheroids. In the middle panel 
of Fig.l, we plot the isodensity curves for this density function: 

P(r) = Pb (l + ~ 2 ) , (2) 

where pb and rt, determine the central density and scale-length of the 
spheroid. 

Bars are non-axisymmetric features with high ellipticities. We will use 
three different bar models. In the first one the bar potential is described by 
a Ferrers ellipsoid (Ferrers [18]) whose density distribution is: 

where m 2 = x 2 /a 2 + y 2 /b 2 . The values of a and 6 determine the shape of 
the bar, a being the length of the semi-major axis, which is placed along 
the x coordinate axis, and b being the length of the semi-minor axis. The 
parameter n measures the degree of concentration of the bar and po rep- 
resents the bar central density. In the right panel of Fig. 1, we plot the 
density function along the semi-major and semi-minor axes of the Ferrers 
ellipsoid with index n = 2, and principal axes a = 6 and b = 1.5. 

We also use two ad-hoc potentials, namely a Dehnen's bar type (Dehnen 
[19]) and a Barbanis-Woltjer (BW) bar type (Barbanis & Woltjer [20]) to 
compare to the results obtained with the Ferrers ellipsoid. The Dehnen's 
bar potential has the following expression: 



$i(r,0) = -ie Wo 2 cos(2#i < 



( ( r\ n 
2-(-) ,r<« 

fa\ n 

(-) . r>a, 



where the parameter a is a characteristic length scale and vo is a character- 
istic circular velocity. The parameter e is related to the bar strength. The 
BW potential has the expression: 

$ 2 (r, 9) = eVr(n - r) ao&{29), (5) 
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where the parameter n is a characteristic scale length and e is related to 
the bar strength. 




Fig. 1. Characteristics of the components. Left panel: Density function of the Kuzmin- 
Toomre disc (red solid line) with = 0.75 and Vd = 1-5. Middle panel: Isodensity curves 
for the spherical distribution representing the spheroid with parameters = 0.3326 and 
Pb = 23552.37. Right panel: Density along the semi-major axis (black solid line) and 
the semi-minor axis (red dashed line) of a Ferrers bar with n = 2, a = 6, b = 1.5 and 
po = 0.0193. 

The bar-like component rotates anti-clockwise with angular velocity 
Jl p = Q p z, where Q p is a constant pattern speed a . 

3. Equations of motion and dynamics around L\ and L2 

The equations of motion in a frame rotating with angular speed S7 p in 
vector form are 

f = -V* - 2{fl p x f) - fJ p x (fJ p x r), (6) 

where the terms 2S~i p x r and — f2 p x (fl p x r) represent the Coriolis and 
the centrifugal forces, respectively, $ is the potential and r is the position 
vector. We define an effective potential <f> cff = <f> — ^fl 2 (x 2 + y 2 ), then Eq. 
(6) becomes f = — V$ ff — 2(O p x f ), and the Jacobi constant is 

£./ = ^|i-| 2 + <I>e ff , (7) 

which, being constant in time, can be considered as the energy in the ro- 
tating frame. The surface <& e ff = Ej (Ej defined as in Eq. (7)) is called the 
zero velocity surface, and its intersection with the z — plane gives the 
zero velocity curve. All regions in which $ e g > Ej are forbidden to a star 
with this energy, and are thus called forbidden regions. 

For our calculations we place ourselves in a frame of reference corotating 
with the bar, and the bar semi-major axis is located along the x axis. In this 



a Bold letters denote vector notation. The vector z is a unit vector. 



February 2, 2008 15:46 



WSPC - Proceedings Trim Size: 9in x 6in romerogomez 



5 

rotating frame we have five equilibrium points, which, due to the similarity 
with the Restricted Three Body Problem, are also called Lagrangian points 
(see left panel of Fig. 2). The points located symmetrically along the x 
axis, namely L\ and L 2 , are linearly unstable. The ones located on the 
origin of coordinates, namely L3, and along the y axis, namely L4 and L$, 
are linearly stable. The zero velocity curve defines two different regions, 
namely, an exterior region and an interior one that contains the bar. The 
interior and exterior regions are connected via the equilibrium points (see 
middle panel of Fig. 2). Around the equilibrium points there exist families 
of periodic orbits, e.g. around the central equilibrium point the well-known 
x\ family of periodic orbits that is responsible for the bar structure. 

The dynamics around the unstable equilibrium points is described in 
detail in Romero-Gomez et al. [21]; here we give only a brief summary. 
Around each unstable equilibrium point there exists a family of periodic 
orbits, known as the family of Lyapunov orbits (Lyapunov [22]). For a given 
energy level, two stable and two unstable sets of asymptotic orbits emanate 
from the corresponding periodic orbit, and they are known as the stable and 
the unstable invariant manifolds, respectively. The stable invariant manifold 
is the set of orbits that tends to the periodic orbit asymptotically. In the 
same way, the unstable invariant manifold is the set of orbits that departs 
asymptotically from the periodic orbit (i.e. orbits that tend to the Lyapunov 
orbits when the time tends to minus infinity), as seen in the right panel of 
Fig. 2. Since the invariant manifolds extend well beyond the neighbourhood 
of the equilibrium points, they can be responsible for global structures. 

In Romero-Gomez et al. [23] , we give a detailed description of the role 
invariant manifolds play in global structures and, in particular, in the trans- 
fer of matter. Simply speaking, the transfer of matter is characterised by 
the presence of homoclinic, heteroclinic, and transit orbits. 

Homoclinic orbits correspond to asymptotic trajectories that depart 
from the unstable Lyapunov periodic orbit 7 around Li and return asymp- 
totically to it (see Fig. 3a). Heteroclinic orbits are asymptotic trajectories 
that depart from the periodic orbit 7 around Li and asymptotically ap- 
proach the corresponding Lyapunov periodic orbit with the same energy 
around the Lagrangian point at the opposite end of the bar Lj, i ^ j (see 
Fig. 3b). There also exist trajectories that spiral out from the region of the 
unstable periodic orbit, and we refer to them as transit orbits (see Fig. 3c). 
These three types of orbits are chaotic orbits since they fill part of the 
chaotic sea when we plot the Poincare surface of section (e.g. the section 
(x, x) near Li). 
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Fig. 2. Dynamics around the L\ and L2 equilibrium points. Left panel: Position of 
the equilibrium points and outline of the bar. Middle panel: Zero velocity curves and 
Lyapunov periodic orbits around L\ and L^- Right panel: Unstable (in red) and stable 
(in green) invariant manifolds associated to the periodic orbit around L\. In grey, we 
plot the forbidden region. From Romero-Gomez et al. 2006, Astronomy & Astrophysics, 
453, 39, EDP Sciences. 




XXX 



Fig. 3. Homoclinic (a), heteroclinic (b) and transit (c) orbits (black thick lines) in the 
configuration space. In red lines, we plot the unstable invariant manifolds associated to 
the periodic orbits, while in green we plot the corresponding stable invariant manifolds. 
In dashed lines, we give the outline of the bar and, in (b) and (c), we plot the zero 
velocity curves in dot-dashed lines. From Romero-Gomez et al. 2007, Astronomy and 
Astrophysics, 472, 63, EDP Sciences. 

4. Results 

Here we describe the main results obtained when we vary the parameters 
of the models introduced in Sec. 2. One of our goals is to check separately 
the influence of each of the main free parameters. In order to do so, we 
make families of models in which only one of the free parameters is varied, 
while the others are kept fixed. Our results show that only the bar pattern 
speed and the bar strength have an influence on the shape of the invariant 



February 2, 2008 15:46 



WSPC - Proceedings Trim Size: 9in x 6in romerogomez 



7 

manifolds, and thus, on the morphology of the galaxy (Romero-Gomez et 
al. [23]). 

Our results also show that the morphologies obtained do not depend 
on the type of bar potential we use, but on the presence of homoclinic or 
hctcroclinic orbits. If heteroclinic orbits exist, then the ring of the galaxy is 
classified as rRi (see Fig. 4a). The inner branches of the invariant manifolds 
associated to 71 and 72 outline an inner ring that encircles the bar and is 
elongated along it. The outer branches of the same invariant manifolds form 
an outer ring whose principal axis is perpendicular to the bar major axis. If 
the model does not have either heteroclinic or homoclinic orbits and only 
transit orbits are present, the barred galaxy will present two spiral arms 
emanating from the ends of the bar. The outer branches of the unstable 
invariant manifolds will spiral out from the ends of the bar and they will 
not return to its vicinity (see Fig. 4d). If the outer branches of the unstable 
invariant manifolds intersect in configuration space with each other, then 
they form the characteristic shape of R2 rings (see Fig. 4b). That is, the 
trajectories outline an outer ring whose principal axis is parallel to the bar 
major axis. The last possibility is if only homoclinic orbits exist. In this 
case, the inner branches of the invariant manifolds for an inner ring, while 
the outer branches outline both types of outer rings, thus the barred galaxy 
presents an i?ii? 2 ring morphology (see Fig. 4c). 



a - rRl ring structure b - rR2 ring structure c - rFi!R2 ring structure d - Barred spiral structure 




-10 -5 5 10 -10 -5 5 10 -20 -10 10 20 -20 -10 10 20 



Fig. 4. Rings and spiral arms structures. We plot the invariant manifolds for different 
models, (a) rR\ ring structure, (b) ri?2 ring structure, (c) R1R2 ring structure, (d) 
Barred spiral galaxy. From Romero-Gomez et al. 2007, Astronomy and Astrophysics, 
472, 63, EDP Sciences. 



5. Summary 

To summarise, our results show that invariant manifolds describe well the 
loci of the different types of rings and spiral arms. They are formed by 
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a bundle of trajectories linked to the unstable regions around the L\jLi 
equilibrium points. The study of the influence of one model parameter on 
the shape of the invariant manifolds in the outer parts of the galaxy reveals 
that only the pattern speed and the bar strength affect the galaxy mor- 
phology. The study also shows that all the different ring types and spirals 
can be obtained when we vary the model parameters. 

We have compared our results with some observational data. Regarding 
the photometry, the density profiles across radial cuts in rings and spiral 
arms agree with the ones obtained from observations. The velocities along 
the ring also show that these are only a small perturbation of the circular 
velocity. 
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